Quantum phase transitions in alternating transverse Ising chains: 
Rigorous analytical and exact numerical results 
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Abstract 

We examine the ground state properties of the s = | transverse Ising chain with regularly alternating bonds 
and fields using exact analytical results and exact numerical data for long (up to N = 900) and short (N = 20) 
chains. For a given period of alternation the system may exhibit a series of quantum phase transitions, where the 
number of transitions depends on the concrete set of the parameters of the Hamiltonian. The critical behaviour 
for the nonuniform and the uniform chain is the same. 

PACS number(s): 75.10.-b 
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Currently the investigation of quantum phase transitions is a very active field of physical research jl], |j . Among 
others the quantum magnets represent an interesting class of systems studied theoretically and experimentally 
where quantum phase transitions appear || . The one-dimensional s = \ Ising model in a transverse field (the 
transverse Ising chain) described by the Hamiltonian 

ff = X>n+£ 2/ «+l (1) 

n n 

( s n — \ a ni a n are the Pauli matrices attached to the site n, f2 is the on-site transverse field, 21 is the exchange 
interaction between neighbouring spins) is known to be the simplest system exhibiting a quantum phase transition. 

For positive the system exhibits a zero temperature transition from the ordered (quantum Ising) phase (for 
ft < Q c ) to the disordered (quantum paramagnetic) phase (for fl > fi c ) at a critical point Q c — \I\ @. (Notice, 
that the same kind of transition takes place at £l c — — 1/|.) There has been a great deal of theoretical work on the 
ground state (GS) properties of the transverse Ising chain ([!]). The knowledge of the eigenvalues and eigenfunctions 



of the Hamiltonian (|l|) |3|, |2j makes the problem accessible for exact analysis. Let us briefly summarize some 
of the main exact results for the GS and the excitation gap. (i) At the critical transverse field f2 c the gap in the 
energy spectrum A closes as A ~ \Sl — f2 c |. (ii) The GS correlation function (sjs| +n ) tends to m x2 ^ as n — > oo 

for < ft < il c with the longitudinal (Ising) magnetization m x = \ {\ — (jr^J ^ which plays the role of the 
order parameter, (iii) For O > O c the correlation function (sjsj +n ) decays exponentially to zero for large n. (iv) 
Close to the critical point tt c the transverse magnetization m z and the static transverse susceptibility \ z behave 
as m z ~ (f2 — fl c ) In |fi — fi c | and x z ~ hi |f2 — fl c \, respectively, (v) At the critical point = £l c the correlation 
length diverges and (sjsj +n ) becomes scale invariant, i.e., it goes to zero for large n according to a power-law 
((s|sj +n ) ~ n~i). (vi) As can be seen from the above results the set of the critical exponents for the quantum 
phase transition in dimension one is identical to that for the thermal phase transition of the square-lattice Ising 
model. 

Recently several authors (see, e.g., Refs. ]?], [|, pi [r], [HJ [l2|, [r], |l4|]) have addressed the question, how deviations 
from the pure uniform crystalline system (e.g., disorder, aperiodically or regularly varying parameters) influence 
the properties of quantum spin systems. In what follows we study this question for the transverse Ising model with 
fields and interactions which vary regularly from site to site with a finite period p, i.e., we change in (|l|) Q — > Q n , 
I — > 7„ and the sequence of parameters along the chain is IR/if^-^ ■ ■ • £l p I p £l lli^k ■ ■ • Q P I P ■ ■ ■ ■ In particular, 
we examine the influence of the so-defined regular inhomogeneity on the existence and the nature of the quantum 
phase transition. In order to study the properties of the regularly alternating transverse Ising chain we use i) exact 
analytical results for thermodynamic quantities obtained by means of continued fractions [ fl5[ , ii) exact analytical 
results for the GS wave function for special parameter values, iii) exact numerical data for long chains of up to 
N = 900 spins |r| , and iv) exact diagonalization data for short chains of N — 20 spins. 

As described in Jl5| the thermodynamic quantities of the regularly alternating transverse Ising chain can be 
calculated rigorously for any finite period of inhomogeneity p. After using the Jordan- Wigner fcrmionization the 
transverse Ising chain is described by noninteracting spinless fermions with the Hamiltonian H = ^2 k A/. (i^Vk ~ 5) 
and the distribution of the square of energy of the elementary excitations R(E 2 ) = J^k $ (E 2 ~ Aj?) can be found 
exactly by means of continued fractions. The resulting density of states has the form 

f -^ |z ;- l(g2)l , if .4 2p (£ 2 )>o, 

I 0, otherwise, 

where Z p -i(E 2 ) and A2 P (E 2 ) = — rTj=i (E 2 — are polynomials of (p — l)th and (2p)th orders, respectively, 
and the a,- > are the roots of A 2p (E 2 ). For example, Z X {E 2 ) = 2E 2 - Of - Q| - if - 1%, Ai(E 2 ) = m\VL 2 2 lll% - 



(E 4 - (Of + n\ + I 2 + 7|) E 2 + n\VLl + (for details of calculation see (15)). R(E 2 ) permits to obtain the 

energy gap A for the spin model, since the smallest root of A2 P {E 2 ) yields the squared minimal excitation energy 
of the spin system. Knowing R{E 2 ) we immediately get all GS quantities which follow from the GS energy per site 
eo = — J °° dEE 2 R(E 2 ). Assuming that fi n = fi + A51„ we obtain the GS transverse magnetization m z = ^| and 
the GS static transverse susceptibility \ z — The thermodynamic quantities at nonzero temperature follow 

from the Helmholtz free energy per site / = —2kT / °° dEER(E 2 ) In (2 cosh jf^)- 

Within the analytical approach based on continued fractions it is not possible to calculate the spin correlation 
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functions. However, we can calculate the GS spin correlation functions (s"s" +n ) (a = x, z) numerically for finite 
chains. It is a specific feature of the considered quantum spin model that the 2 N x 2 N problem inherent in 
the exact diagonalization scheme reduces to a N x N problem || || and therefore numerical exact data for 
rather long chains can be obtained. (Notice, that typical chain length accessible for Lanczos diagonalization of 
quantum Heisenberg systems with low symmetry is about N = 24.) The on-site magnetizations can be calculated 
from the correlation functions using the relation Tn^m^ = lim r ^ 00 (s" i s" 2+) . p ), where 1 <C ji <C N is one of 
the p consecutive numbers sufficiently far from the ends of chain and j% — ji — 0, 1, ... ,p — 1. Assuming that 
( s f s ?+rp) ~ ( s< j) ( s< j+r P ) ~ ex P ( — f«) f° r l ar S e r (1 <ti j,j +rp -c TV) we can determine the correlation length £ a 
from the exact numerical data. 

Let us discuss the effects of regular alternation on the existence, the position and the critical properties of the 
quantum phase transition in the transverse Ising chain in more detail. For a certain set of parameters we identify 
the points of quantum phase transition by looking for vanishing excitation gap A. As mentioned above A is 
determined by the smallest root of A2 P (E 2 ) and therefore the location of a quantum phase transition is determined 
by A2 P (0) = 0. For period p = 2 the corresponding condition ^(0) = yields 

ftifi 2 - ±IiJ»- (3) 

It should be mentioned that a condition for closing of the excitation gap for an inhomogeneous transverse Ising 
chain was already reported in Ref. |Tij| . Eq. (0) agrees with that result for the regularly alternating chain of 
period 2. In what follows we consider as an example a chain with Sl± % = O ± AO, AO > 0. We find either two 



characteristic fields O c = ±i/AO + I/1I2I if AO < y/\Iih\ or four characteristic fields O c = ±y AO ± \Iih\ if 



AO > y^/iT^I p^q| . This is illustrated in Figs, la, Id where the gap A in dependence on (dotted curves) for 
the chains with |7i| = |7 2 | = 1 and Afl — 0.5 (Fig. la) and Ail — 1.5 (Fig. Id) is shown. Thus, in case of small 
strength of transverse-field inhomogeneity Afl only quantitative deviations from the homogeneous case may be 
expected. Contrary, strong inhomogeneity AO increases the number of quantum phase transitions |H| . Using the 
explicit expression for Aa(E 2 ) we find that the smallest root a,j decays as ~ (fi — il c ) 2 for — > fl c . Therefore, the 
energy gap A closes as ~ |0 — fi c | for — > O c , whereas the GS energy eo contains the term ~ (0 — O c ) 2 In |0 — O c | 
and hence m z and \ z contain the terms ~ (Q — Q c ) In \fl — f2 c | and ~ In |fi — f2 c |, respectively. This means that 
the critical behaviour of the regularly alternating transverse Ising chain is identical to that of the homogeneous 
chain. In Figs, la, Id we show m z vs. (solid curves) and \ z vs. fl (dashed curves) for Afl = 0.5 (Fig. la) and 
AO = 1.5 (Fig. Id). We emphasize, that the quantum phase transition is not suppressed by the deviation of the 
uniformity in kind of a regular alternation of bonds and fields. Even the appearance of additional transitions due 
to field alternation is possible. 

In addition to the exact results for thermodynamic quantities we can obtain exact expressions for the GS 
|GS) of the considered chain at special values of O. Obviously, we have |GS) = FT | Tn) (|GS) = Y[ n I In)) f° r 
-oo (O - +00), |GS) = . . . I U| Tn+i) ... for O = and AO » \h\, \I 2 \ and |GS) = U n 75 (I in) + I Tn» 
or |GS) = Yin 772 (I i") ~ I in)) ^ or ^ = = ^ ( we nave assumed the ferromagnetic sign of the exchange 
interactions). Nontrivial but solvable is the case O = =fAO, i.e., Oi = 0, O2 = — 2AO or Oi = 2AO, O2 = 0. 
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Supposing > the GS for fi n = 0, fi n +i = — 2A17 is given by 

|GS> = U+R + ...( Cl \ i n+1 )\ 4„ +2 >+c 2 | T„+i)| Tn+a»(ci| k+3>l k+ 4 )+c 2 | T„+3>| T„ +4 » • • • , (4) 

U + = . . . exp (-ms x n s y n+1 ) exp (-iTrs^^) 



R + = ... exp (-i - a* J exp (-i - s*^ 

1 I+^P + Afi 2 1 Aft 7i 

c i = — ^= , . c 2 



\/2 V/ 2 + /V/ 2 + Aft 2 + Aft 2 ' V2 V/ 2 + V/ 2 + Aft 2 + Aft 2 ' 2 

To obtain (Q) we note that the transformed Hamiltonian RUHU + R + for AT — > oo (when boundary terms become 
negligible) describes a system of noninteracting two-site clusters each with the Hamiltonian I\sf l+1 + / 2 sf l+2 — 
4Afts£ +1 s^ +2 . For ft n = 2Aft, ft„+i = the GS is again given by (Q) with the change n — ► n — 1, Aft — > —Ail. 
Knowing |GS) we can calculate the xx and zz spin correlation functions. For example, for ft = —Aft according to 
Eq. (j|) we find 

( s n+l s n+2) — ( S n+l S n+4) = ■ ■ • = { s n+2 s n+'il = ( s n+2 s n+h) = ■ ' • = 7 ( C 2 



/ 2 2\ 2 /ar xv/ararv ^ 
— lyC 2 — CjJ , {S n+2 S n+ 4) = \S n+ 2 s n+6) = ' ' ' = 4' 

( S n+l S n+2) = ( S n+l S n+4) = • •• = ( S ri+2 S n+3) = ( S n+2 S n+5) = ■ ■ ■ = { S 7i+2 S 7i+4> = ( S n+2 S n+6) = ■ • • = 0, 



s n+l*n+3) _ (*n+l*n+5) • • • "7 c l) ! - • ~ •."...-•.'"•••Hi.' 



'n+l^n+3/ — \°n+l''n+5 



> = ... = (c lC2 ) 2 . (5) 



For ft = Aft the correlation functions are given by (||) with the change n — > n — 1. 

Now we discuss the exact numerical results for finite chains presented in Fig. 1 and Fig. 2. The order parameters 



(longitudinal sublattice magnetizations) \m x \ versus transverse field ft are shown in Figs. 1 and 2 for Aft < y / ]ii7 2 | 



(Figs, lb, 2a (open symbols)) and Ail > \J\I1I2\ (Figs, le, 2a (full symbols)). In accordance with the analytical 
results for A, m z and \ z the numerical data show the existence of either two phases (quantum Ising phase for 



|Q| < J Aft + |ii/ 2 | and strong-field quantum paramagnetic phase otherwise) or three phases (low-field quantum 
paramagnetic phase for |ft| < J Ail 2 — {hhl, quantum Ising phase for J Aft 2 — \I\h\ < jO| < \J Aft 2 + |/i/ 2 | 
and strong-field quantum paramagnetic phase otherwise) . The longitudinal on-site magnetizations m x are nonzero 
in quantum Ising phases and become strictly zero in quantum paramagnetic phases. The transverse magnetization 
to 2 in quantum paramagnetic phases is almost constant being in the vicinity of zero in the low-field phase and of 
saturation value in the strong-field phase thus producing plateau-like steps in the dependence m z versus ft. The 
correlation length £ x (Figs, lc, If, 2b (full symbols)) illustrates that the transition between different phases is 
accompanied by the divergency of £ x . As can be seen in Figs, lb, le short chains are sufficient to reproduce the 
correct order-parameter behaviour in the quantum Ising phase away from the quantum critical point, but not in the 
quantum paramagnetic phases. The chain length of up to A^ = 900 sites is clearly sufficient to see a sharp transition 
in the order parameter (Figs, lb, le) and even to extract reliable results for £ x from the long-distance behaviour 
of (s x s x +r ) as can be seen from the data presented in Figs, lc, If and in Fig. 2b. We note that the numerical data 
for \m x \ and £ x at ft = =FAft coincide with the exact expression (||). In particular, the values of the longitudinal 
sublattice magnetizations for these fields are | and \ \c 2 — c\\. It should be remarked that inhomogeneity inherent 
in the model is regular but not random and there is no frustration; as a result we do not observe any indication of 
a glassy behaviour. 
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Finally, the low-temperature dependence of the specific heat c at different (> 0) obtained from the exact 
analytical expression for the free energy (see above) confirms the existence of either two phase transitions (Fig. 
3a) or four phase transitions (Fig. 3b) depending on the relationship between Afi and \J\I1I2\. At the quantum 
phase transition points the spin chain becomes gapless and c depends linearly on T. The ridges seen in Figs. 3a, 
3b (which correspond to maxima in the dependence c vs. f2 as T varies) single out the boundaries of quantum 
critical regions O . These boundaries correspond to a relation A ~ kT that can be checked by comparison with the 
data for A vs. O reported in Figs, la, Id. As can be seen from Fig. 3 the c(T) behaviour for O slightly above or 
below VL C changes crossing the boundaries of quantum critical regions. Furthermore we notice that for A£7 = 1.5 
an additional low-temperature peak appears in c(T). 

To end up, let us turn to a chain of period 3 for which the analytical and numerical calculations presented above 
can be repeated. The quantum phase transition points follow from the condition Aq(0) — 0, i.e., f2i 172^3 = ±ii/2 -fa- 
Depending on the parameters either two, four, or six quantum phase transitions are possible. This behaviour is 
illustrated in Fig. 4 for a chain of period 3 with |ii| = \I 2 \ = \I 3 \ = 1, fi n = f2 + AQ„, Afii + Ail 2 + Ail 3 = 0. 

To summarize, we present exact results for the thermodynamics of a regularly alternating s = | Ising chain 
in a transverse field. Furthermore, for special parameter values also the exact GS and spin correlation functions 
are given. For parameters for which the correlation functions are not accessible for rigorous analysis we present 
exact numerical results for finite but large systems. From these numerical data we can extract the order parameter 
and the correlation length in high precision. We find, that the quantum phase transition present in the uniform 
chain is not suppressed by deviation from uniformity in kind of a regular alternation of bonds and fields. On the 
contrary, field alternation may lead to the appearance of additional transitions. The number of phase transitions 
for a given period of alternation strongly depends on the precise values of the parameters of the model. For the gap 
A, the GS energy eo, the transverse magnetization m z and the static transverse susceptibility \z we can examine 
the critical behaviour exactly and find the same critical indices for the nonuniform and the uniform chain. The 
presented study can be extended for examining the ground state properties of s = 5 transverse Ising model on 
one-dimensional superlattices. 

The present study was partly supported by the DFG (projects 436 UKR 17/7/01 and 436 UKR 17/7/02). O. 
D. acknowledges the kind hospitality of the Magdeburg University in the summer of 2002. O. D. and T. K. were 
supported by the STCU under the project #1673. 
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FIGURE CAPTIONS 

Figure 1. The GS properties of the transverse Ising chain of period 2 with |ii| = |J 2 | = 1 and H12 = fi ± Ail 
in dependence on the transverse field Q: Energy gap A (dotted curves) (a,d), transverse magnetization m z (solid 
curves) (a,d), static transverse susceptibility \ z (dashed curves) (a,d), longitudinal sublattice magnetizations |mj| 
(triangles and solid curves correspond to the data for N = 20 and N = 600, respectively) (b,e), and correlation 
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length £ x (triangles and squares correspond to the data for N = 300 and TV = 600, respectively) (c,f) for AO = 0.5 
(a,b,c) and AO = 1.5 (d,e,f). 

Figure 2. The GS properties of the transverse Ising chain of period 2 with |ii| = \I 2 \ = 1 and 0^ 2 = fi ± AO 
in dependence on the transverse field O: Longitudinal sublattice magnetizations |mj| (a) and correlation length 
£ x (b) for chain length of N = 300 (triangles), N = 600 (squares), N = 900 (circles) sites and AO = 0.99 (open 
symbols), AO = 1.01 (full symbols). 

Figure 3. The low- temperature dependence of the specific heat for the transverse Ising chain of period 2 with 
l^i I = l-k | = 1 and fi 1)2 = fi± AO; AO = 0.5 (a), AO = 1.5 (b). 

Figure 4. Number of quantum phase transitions present in the transverse Ising chain of period 3 with = 
1*2 1 = |*3 1 = 1 and 0„ = O + AO„, AOi + A0 2 + A0 3 = 0. The dark, grey, or light regions correspond to 
parameters where the system exhibits two, four, or six quantum phase transitions, respectively. 
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